AD-A241  015 


TECHNICAL  REPORT  BRL-TR-3265 


BRL 


U  ( 


t  * 


C  i  ■ 


$1 

rr-.  vV**' 


TRANSONIC  NAVIER-STOKES  COMPUTATIONS 
FOR  A  SPINNING  BODY  OF  REVOLUTION 


JUBARAJ  SAHU 


SEPTEMBER  1991 


<£> 

I 

M 

U1 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


U.S.  ARMY  LABORATORY  COMMAND 


BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MARYLAND 


NOTICES 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the  originator. 

Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical  Information  Service, 
U.S.  Department  of  Commerce,  5285  Port  Royal  Road,  Springfield,  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position, 
unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers’  names  in  this  report  does  not  constitute  indorsement 
of  any  commercial  product. 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OM8  No  0704-0188 


Puolic  reoortir.g  ouraen  tor  thii  collection  ot  ntorirMtion  s  e^timaieo  to  i«er*qe  •  ■'our  oer  'esoorse.  nciuamg  the  time  to r  reviewing  instruction*.  searching  ?«  sting  a<*ta  sources, 
gathering  jnd  mamtjtnmg  the  data  needed,  and  comoietmg  ana  reviewing  th*  .oiieaion  ot  information  Send  comments  regarding  this  Durden  estimate  or  any  ither  asoect  ot  trus 
collection  ot  information,  nciudmg  suggestions  tor  reducing  this  ouroen  to  Washington  neadauarters  Services.  Directorate  for  information  Ooerations  and  *eocm.  ‘  2 1 5  .e'terscn 
Davis  Highway.  Su.te  *204.  Arlington,  *  a  222Q2-J302  and  to  the  Oftice  o*  Management  and  Budget.  Pioerwor*  Reduction  Project  (07C4-Q188).  Wasnmgton.  DC  20S03 


1.  AGENCY  USE  ONLY  (Leave  blank) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  OATE 

September  1991 


3.  REPORT  TYPE  AND  DATES  COVERED 
Final,  July  1989  -  July  1991 


S.  FUNDING  NUMBERS 


TRANSONIC  NAVIER- STOKES  COMPUTATIONS  FOR  A  SPINNING  BODY  OF 
REVOLUTION 


6.  AUTHOR(S) 


JUBARAJ  SAHU 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 


1L161102AH43 
61 102A-00-001  AJ 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME{S)  AND  ADDRESS(ES) 

US  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  Maryland  21005-5066 


10.  SPONSORING.  MONITORING 
AGENCY  REPORT  NUMBER 


BRL— TR-3265 


12a.  DISTRIBUTION/ AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13.  abstract  (Maximum  200 words)  ^  zonal,  implicit,  time-marching  Navier-Stokes  computational 
technique  has  been  used  to  compute  three  dimensional  transonic  flow  fields  over  a 
projectile.  Flow  field  computations  have  been  performed  at  M  =  0.94  for  spin  rates 
of  0  and  4900  rpm  and  at  angles  of  attack,  alpha  »  0,  4,  and  10  degrees.  All  the 
computations  have  been  performed  on  the  Cray-2  supercomputer.  Details  of  the  flow 
field  such  as  Mach  number  contours  and  surface  pressure  distributions  are  presented. 
Computed  surface  pressures  are  compared  with  available  experimental  data  for  the  same 
conditions  and  the  same  configuration.  Computed  results  show  the  large  circumferen¬ 
tial  pressure  distribution  over  the  boattail  region  as  well  as  the  nonlinear  effect 
of  the  angle  of  attack.  Aerodynamic  force  and  moment  coefficients  (normal  force, 
pitching  moment,  Magnus  force,  and  Magnus  moment)  have  been  obtained  from  the 
computed  pressures  and  are  compared  with  the  data.  The  computed  results  are 
generally  in  good  agreement  with  the  data  at  low  angle  of  attack  for  both  nonspinning 
and  spinning  conditions.  At  high  angle  of  attack  the  agreement  is  good  for  the 
nonspinning  case  and  is  less  satisfactory  for  the  spinning  case. 


J4.  SUBJECT  terms 
Aerodynamics 

Three  Dimensional  Flow 

Spinning 

Transonic  Flow 


17  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


■jSM  75*10-0’  -280-5500 


Magnus  Effect 
Angle  of  Attack 
Computations 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


15.  NUMBER  OF  PAGES 

37 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


^tardara  -■ 


:98  Rev  :-89) 


Intentionally  left  blank. 


ii 


TABLE  OF  CONTENTS 


Page 


LIST  OF  FIGURES .  v 

LIST  OF  TABLES  .  vi 

1.  INTRODUCTION .  1 

2.  GOVERNING  EQUATIONS  AND  SOLUTION  TECHNIQUE .  2 

2.1  Governing  Equations .  2 

2.2  Numerical  Technique .  4 

2.3  Composite  Grid  Scheme .  5 

3.  MODEL  GEOMETRY  AND  EXPERIMENT .  6 

4.  COMPUTATIONAL  GRID  AND  COMPUTER  RESOURCE .  7 

5.  RESULTS  .  8 

6.  CONCLUDING  REMARKS .  27 

7.  REFERENCES  .  29 

LIST  OF  SYMBOLS .  31 

DISTRIBUTION  LIST .  33 


A  .  it  It  B  i  tf 

;  a”'-'-’  <JfU4I 

v  t :  t»jb 

,  J-lSLi/lc&l  loc. 


■ .  *'  rt  bt> t.  i*n/ 

i  iity  Cc4«b 

j  iAitm'I  h  on  /  • 

'  ^  ■•f*  i  'i  i 


!f\A  ! 


iii 


INTENTIONALLY  LEFT  BLANK. 


i 


IV 


Figure 


LIST  OF  FIGURES 


Page 

1  Model  geometry  .  6 

2  Computational  grid  (wind  and  lee  side) .  8 

3  Expanded  view  of  the  base  region  grid .  9 

4  Mach  number  contours,  Moo=0.94,  a  =  0.0° .  10 

5  Mach  number  contours,  Moo=0.94,  a  =  4.0° .  11 

6  Mach  number  contours,  Moo=0.94,  a  =  10.0° .  11 

7  Longitudinal  surface  pressure  distribution,  Moo=0.94,  a  =  0.0° .  12 

8  Longitudinal  surface  pressure  distribution,  Moo=0.94,  a  =  4.0° .  12 

9  Longitudinal  surface  pressure  distribution,  Moo=0.94,  a  =  10.0° .  13 

10  Circumferential  surface  pressure  distribution,  P  =  0,  Moo=0.94,  X/D=4.76  .  14 

11  Circumferential  surface  pressure  distribution,  P  =  0,  Moo=0.94,  X/D=5.11  .  14 

12  Circumferential  surface  pressure  distribution,  P  =  0,  Moo=0.94,  X/D=5.30  .  15 

13  Circumferential  surface  pressure  distribution,  P  —  0,  Moo=0.94,  X/D=5.43  .  15 

14  Circumferential  surface  pressure  distribution,  P  —  0,  Moo=0.94,  X/D=5.50  .  16 

15  Circumferential  surface  pressure  distribution,  P  =  4900  rpm,  Moo =0.94, 

X/D=4.76 .  17 

16  Circumferential  surface  pressure  distribution,  P  =  4900  rpm,  Moo=0.94, 

X/D=5.11 .  17 

17  Circumferential  surface  pressure  distribution,  P  =  4900  rpm,  Mcc,=0.94, 

X/D=5.30 .  18 

18  Circumferential  surface  pressure  distribution,  P  =  4900  rpm,  Moo=0.94, 

X/D=5.43 .  18 

19  Circumferential  surface  pressure  distribution,  P  =  4900  rpm,  Moo=0.94, 

X/D=5.50 .  19 

20  Circumferential  surface  pressure  distribution,  Moo=0.94,  a  =  4.0°,  X/D=4.76  20 

21  Circumferential  surface  pressure  distribution,  Moo=0.94,  a  =  10.0°,  X/D=4.76  20 

22  Circumferential  surface  pressure  distribution,  Mco= 0.94,  a  =  4.0°,  X/D=5.43  21 

23  Circumferential  surface  pressure  distribution,  Moo=0.94,  a  =  10.0°,  X/D=5.43  21 

24  Normal  force  coefficient  vs.  angle  of  attack,  Mcc=0.94 .  22 


V 


25  Pitching  moment  coefficient  vs.  angle  of  attack,  M0 <>=0.94 .  23 

26  Development  of  Magnus  force  along  the  projectile,  M(X= 0.94 .  24 

27  Magnus  moment  coefficient,  Moo  =0.94 . 25 


28  Circumferential  surface  pressure  distribution,  Moo=0.94,  a  —  10.0°,  X/D=4.76  26 

29  Circumferential  surface  pressure  distribution,  Moo=0.94,  a  =  10.0°,  X/D=5.50  27 


LIST  OF  TABLES 


Table  Page 

1  Analysis  of  Nonlinear  Magnus  moment .  25 


VI 


1.  INTRODUCTION 


A  spinning  projectile  in  flight  at  angle  of  attack  produces  the  Magnus  effect.  Many 
Army  projectiles  are  slender,  spin-stabilized  bodies  of  revolution  and  at  angle  of  attack  are 
thus,  subject  to  Magnus  forces  and  moments.  These  forces  and  moments  are  generated  by 
the  distortion  or  asymmetry  of  the  boundary  layer  with  respect  to  the  plane  of  angle  of 
attack.  Although  the  Magnus  force  is  only  a  10th  to  a  100th  of  the  normal  force,  it  is  a 
critical  parameter  in  determining  the  dynamic  stability  of  shell.  In  addition,  the  critical 
aerodynamic  behavior  occurs  in  the  transonic  speed  regime,  0.9  <  M  <  1.1,  where  the 
aerodynamic  coefficients  have  been  found  to  increase  by  as  much  as  100  %.  This  flight 
regime  is  both  experimentally  and  computationally  difficult.  Thus,  only  a  limited  number 
of  experimental  and  computational  studies  have  been  made  for  spinning  projectiles  in  this 
speed  regime. 

Numerical  capabilities  based  on  Navier-Stokes  computational  technique  have  been  de¬ 
veloped  and  used  to  compute  three  dimensional  transonic  flowfields  (Nietubicz,  Sturek, 
and  Heavey  1983;  Sahu  1987,  1988;  Sahu  and  Steger  1987;  Sahu  and  Nietubicz  1989). 
Earlier  three  dimensional  computations  for  spinning  (Nietubicz,  Sturek,  and  Heavey  1983) 
and  non-spinning  (Sahu  1987)  projectiles  at  transonic  speed  suffered  from  inadequate  grid 
resolution.  Recently,  improved  three  dimensional  (3D)  predictions  (Sahu  and  Steger  1987; 
Sahu  iy8b;  Sahu  and  Nietubicz  J.9fi9)  for  non-spmning  projectiles  have  been  obtained  using 
the  Cray-2  supercomputer  and  a  more  advanced  upwind  flux-split  algorithm.  The  critical 
aerodynamic  behavior  in  pitching  moment  coefficient  was  successfully  predicted  using  a 
zonal  (Sahu  and  Steger  1987),  implicit,  time-marching  Navier-Stokes  scheme.  This  zonal 
scheme  preserves  the  base  corner  and  allows  better  modeling  of  the  uuec  region  flow.  The 
present  research  is  a  further  extension  and  application  of  this  technique  to  a  spinning 
projectile  at  angles  of  attack  up  to  10  degrees. 

The  basic  configuration  used  in  this  study  has  a  3.0  caliber  nose,  2.1  caliber  cylinder 
and  a  7°,  0.5  caliber  boattail.  The  experimental  model  was  sting  mounted  from  the  base. 
As  mentioned  out  earlier,  a  simple  composite  grid  scheme  has  been  used  for  accurate 
modeling  of  the  base  corner.  Numerical  flowfield  computations  have  been  performed  at 
M  =  0.94  for  spin  rates  of  0  and  4900  rpm  and  at  angles  of  attack,  a  =  0,4,  and  10°  . 
All  the  computations  have  been  performed  on  the  Cray-2  supercomputer.  Details  of  the 
flowfield  such  as  Mach  number  contours  and  surface  pressure  distributions  are  presented. 
Computed  surface  pressures  are  compared  with  available  experimental  data  for  the  same 
conditions  and  the  same  configuration.  Computed  results  show  the  h.gc  change  in  the 
circumferential  pressure  distribution  over  the  boattail  region  as  well  as  the  non-linear  effect 
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of  the  angle  of  attack.  Aerodynamic  force  and  moment  coefficients  (normal  force,  pitching 
moment  coefficient,  Magnus  force,  and  Magnus  moment)  have  been  obtained  from  the 
computed  pressures  and  are  compared  with  the  data. 


2.  GOVERNING  EQUATIONS  AND  SOLUTION  TECHNIQUE 

The  complete  set  of  time-dependent  thin-layer  Navier-Stokes  equations  is  solved  nu¬ 
merically  to  obtain  a  solution  to  this  problem.  The  numerical  technique  used  is  an  implicit 
finite  difference  scheme.  Although  time-dependent  calculations  are  made,  the  transient 
flow  is  not  of  primary  interest  at  the  present  time.  The  steady  flow,  which  is  the  desired 
result,  is  obtained  in  a  time  asymptotic  fashion. 


2.1  Governing  Equations.  The  complete  set  of  three  dimensional,  time  depen¬ 
dent,  generalized  geometry,  thin-layer,  Navier-Stokes  equations  for  general  spatial  coordi¬ 
nates  £,  77,  £  can  be  written  as  (Pulliam  and  Steger  1982): 

dTq  +  fyF  +  dnG  +  dcH  =  Re-'d^S  (1) 

where 

£  *■  '  M) 

n  - 

C  =C  (x,y,z,t) 

T  =  t 

and 


-  longitudinal  coordinate 

-  circumferential  coordinate 

-  nearly  normal  coordinate 

-  time 


p 
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pu 
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.  (e  +p)V  -  rjtP  J 
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(2) 


2 


ana  where 


0 


/XCx  +  Cy  +  C)u<  +  ^-(Cx^C  +  C><  +  Czwc  )Cx 

MCx2  +  Cy  +  Cz)v(  +  3^rU<  +  C sTC  +  Czu,c)Cv 

C  _  1  /XCx  +  Cy  +  C;)^<  +  ^((xU<  +  Cyu<  +  Gtr<  )C- 

J  J 

{(Cx  +  Cy  +  Cx)i^(^2  +  l>2  +  U>2)( 

+ — ^ — ] 

+  ^(Cxu  +  Cyv  +  Cz^XC  xuC  +  C»u<  +  Cx^c)} 

In  equation  (1),  the  thin-layer  approximation  is  used  and  the  viscous  terms  involving 
velocity  gradients  in  both  the  longitudinal  and  circumferential  directions  are  neglected. 
The  viscous  terms  are  retained  in  the  normal  direction,  f,  and  are  collected  into  the  vector 
5.  These  viscous  terms  are  used  everywhere  except  in  the  wake  (or  base  region)  where 
similar  viscous  terms  axe  also  added  in  the  streamwise  direction. 

For  this  computation,  the  diffusion  coefficients  p  and  k  contain  molecular  and  turbu¬ 
lent  parts.  The  turbulent  contributions  are  supplied  through  an  algebraic  eddy-viscosity 
hypothesis  which  has  been  developed  by  Baldwin  and  Lomax  (1978). 

The  velocities  in  the  f,  r/,  and  (  coordinate  directions  can  be  written 

U  =  ft  +  ufx  +  ufy  4-  U’f, 

V  -  Tjt  +  UT)X  4-  rn/y  +  V!TJZ 
W  =  Ct  +  «Cx  +  V(y  +  w(z 

which  represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  retained  as  the  dependent  variables 
and  are  nondimensionalized  with  respect  to  a <*>  (the  free  stream  speed  of  sound).  The  local 
pressure  is  determined  using  the  relation 

p  -  (7  -  l)[e  -  0.5p(u2  -4  v2  +  u;2)]  (4) 
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where  7  is  the  ratio  of  specific  heats.  Density,  p ,  is  referenced  to  p^  and  the  total  energy, 
e,  to  people  ■  The  transport  coefficients  are  also  nondimensionalized  with  respect  to  the 
corresponding  free  stream  variables.  Thus  the  Prandtl  number  which  appears  in  5  is 
defined  as  Pr  =  c^p^/ 

In  differencing  these  equations  it  is  often  advantageous  to  difference  about  a  known 
base  solution  denoted  by  subscript  0  as 


tr(Q-Qo)  +6((F  -  Fo)  +  Sn{G -Go)  +  6^H  -  H0)~  Re~l6c(S  -  So)  (5) 
=  -8tQo  -  d,F0  -  dnG0  ~  d<Ho  +  Re~ldcS0 

w'here  6  indicates  a  general  difference  operator,  and  d  is  the  differential  operator.  If  the 
base  state  can  be  properly  chosen,  the  differenced  quantities  can  have  smaller  and  smoother 
variation  and  therefore  less  of  a  differencing  error  (Pulliam  and  Steger  1982). 


2.2  Numerical  Technique.  The  implicit,  approximately-factored,  scheme  for  the 
thin  layer  Navier  Stokes  equations  using  centred  differencing  in  the  7  and  £  directions  and 
upwinding  in  £  is  written  in  the  following  form 


[/  +  hSl(A+)n  +  h6 c_Cn  -  h.Re~l6(J~l Mn J  -  D, |c 
x  [/  +  h6'(A~)n  +  h6vBn  -  A|„]  AQ"  = 
-A t{8b([(F+r  -  F+]  +  6{l(F-)n  ~  F~]  +  6,(6"  -  G„) 
+6c(Hn  -  tfoo)  -  Re-'SdS*  -  S oo)}  -  De(Qn  -  Q^) 


(6) 


where  h  —  A t  or  (A<)/2  and  the  free  stream  base  solution  is  used.  Here  6  is  typically  a 
three  point  second  order  accurate  central  difference  operator,  6  is  a  midpoint  operator  used 
with  the  viscous  terms,  and  the  operators  6^  and  6^  are  backward  and  forward  three-point 
difference  operators.  The  flux  F  has  been  eigensplit  and  the  matrices  .4,B,C,  and  M 
result  from  local  linearization  of  the  fluxes  about  the  previous  time  level.  Here  J  denotes 
the  Jacobian  of  the  coordinate  transformation.  Dissipation  operators,  De  and  Dl  are  used 
in  the  central  space  differencing  directions. 

The  smoothing  terms  used  in  the  present  study  are  of  the  form: 


De],  =  (A  t)J~'[e26p(B)06  +  e46^6%J 
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D,  I,  =  (A  t)J-l[e2Sp(B)/3S  +  2.oe46p(B)6}\r,J 


where  (3  =  and  where  p(B )  is  the  true  spectral  radius  of  B.  The  idea  here  is  that 

the  fourth  difference  will  be  tuned  down  neax  shocks,  that  is,  as  (3  gets  large  the  weight  on 
the  fourth  difference  drops  down  while  the  second  difference  tunes  up. 

For  simplicity,  most  of  the  boundary  conditions  have  been  imposed  explicitly  (Sahu 
1987).  For  the  non-spinning  cases,  the  wall  temperature  is  specified  on  the  body  surface  and 
the  no-slip  boundary  condition  is  used  at  the  wall.  The  pressure  at  the  wall  is  calculated 
by  solving  a  combined  momentum  equation.  Free  stream  boundary  conditions  are  used  at 
the  inflow  boundary  as  well  as  at  the  outer  boundary.  A  symmetry  boundary  condition 
is  imposed  at  the  circumferential  edges  of  the  grid  while  a  simple  extrapolation  is  used 
at  thi  downstream  boundary.  A  combination  of  symmetry  and  extrapolation  boundary 
condition  is  used  at  the  center  line  (axis).  As  mentioned  in  the  introduction,  the  Magnus 
effect  is  produced  by  the  spin  induced  distortion  of  the  boundary  layer.  This  requires  the 
numerical  calculations  to  be  fully  three  dimensional  since  no  plane  of  symmetry  exists.  For 
the  spinning  cases,  spin  or  the  non-dimensional  angular  velocity  is  provided  through  the 
contravariant,  V  term.  In  addition,  a  periodic  implicit  boundary  condition  is  used  in  the 
circumferential  direction.  The  flowfield  is  initially  set  to  free  stream  conditions  everywhere 
and  then  advanced  in  time  until  a  steady  state  solution  is  obtained. 


2.3  Composite  Grid  Scheme.  In  the  present  work,  a  simple  composite  grid  scheme 
(Sahu  1988)  has  been  used  where  a  large  single  grid  is  split  into  a  number  of  smaller  grids 
so  that  computations  can  be  performed  on  each  of  these  grids  separately.  These  grids  use 
the  available  core  memory  one  grid  at  a  time,  while  the  remaining  grids  are  stored  on  an 
external  disk  storage  device  such  as  the  solid  state  disk  device  (SSD)  of  the  Cray  X-MP/48 
computer.  The  Cray-2  has  a  large  incore  memory  to  fit  the  large  single  grid.  However, 
for  accurate  geometric  modeling  of  complex  projectile  configurations  which  include  blunt 
noses,  shaxp  corners  and  base  cavities,  it  is  also  desirable  to  split  the  large  data  base  into 
a  few  smaller  zones  on  the  Cray-2  as  well. 

The  use  of  a  composite  grid  scheme  requires  special  care  in  storing  and  fetching  the  in¬ 
terface  boundary  data  ,  i.e.,  the  communication  between  the  various  zones.  In  the  present 
scheme,  there  is  a  one  to  one  mapping  of  the  grid  points  at  the  interface  boundaries.  Thus, 
no  interpolations  are  required.  Details  of  the  data  storage,  data  transfer  and  other  perti¬ 
nent  information  such  as  metric  and  differencing  accuracy  at  the  interfaces  can  be  found 
in  the  work  of  Sahu  and  Steger  (1988)  and  Sahu  (1988).  This  scheme  has  been  successfully 
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used  by  Sahu  (1988)  to  compute  three  dimensional  transonic  flow  over  projectiles  at  angle 
of  attack  but  without  spin.  The  computed  results  clearly  showed  the  transonic  critical 
aerodynamic  behavior  of  the  pitching  moment  coefficient  observed  in  free  flight.  The  same 
technique  was  also  appned  to  a  more  complicated  projectile  with  base  cavities  (Sahu  and 
Nietubicz  1989)  and  the  wake  flowfields  were  computed  for  these  different  configurations. 
The  present  work  is  a  further  application  of  this  technique  to  the  spinning  projectile  case. 


3.  MODEL  GEOMETRY  AND  EXPERIMENT 

The  computational  accuracy  of  a  numerical  scheme  can  be  established  through  com¬ 
parisons  with  available  experimental  data.  The  model  used  in  the  experiment  and  in  the 
computational  study  is  shown  in  Figure  1.  It  consists  of  a  3.0  caliber  secant-ogive  nose, 
a  2.0  caliber  cylindrical  section,  and  a  0.5  caliber  7°  boattail.  This  configuration  closely 
resembles  the  M549  155mm  artillery  shell.  The  model  has  a  flat  nose  and  is  sting-mounted 
at  the  aft  end  of  the  projectile.  Both  of  these  are  included  in  the  numerical  simulations 
for  a  direct  comparison. 


Figure  1.  Model  geometry. 

Experimental  measurements  (Miller  and  Molnar  1986)  for  this  model  have  been  made 
at  the  NASA  Ames  14  Foot  Transonic  Wind  Tunnel.  The  model  was  tested  at  angles  of 
attack  of  0,  4,  and  10  degrees  and  spin  rates  of  0  and  4900  rpm.  Testing  was  conducted  at  a 
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Mach  number  of  0.94  and  Reynolds  number  of  4  x  106  per  foot.  Circumferential  pressures 
have  been  measured  at  several  longitudinal  positions  on  the  model,  with  special  emphasis 
on  the  cylindrical  and  boattail  sections.  These  pressures  have  also  been  integrated  to 
obtain  the  aerodynamic  forces  and  moments  which  are  available  for  comparison  with  the 
computed  results. 

4.  COMPUTATIONAL  GRID  AND  COMPUTER  RESOURCE 

The  solution  technique  requires  the  discretization  of  the  entire  flow  region  of  interest 
into  a  suitable  computational  grid.  The  grid  outer  boundary  has  been  placed  at  2.5  body 
lengths  upstream  and  surrounding  the  projectile.  The  downstream  boundary  was  placed 
at  2  body  lengths.  Since  the  calculations  are  in  the  subsonic/transonic  regime  the  compu¬ 
tational  boundaries  must  extend  out  beyond  the  influence  of  the  body.  This  ensures  that 
the  boundary  conditions  specified  in  the  flow  code  are  satisfied. 

Since  transonic  flow  fields  are  complex  due  to  the  presence  of  mixed  subsonic  and 
supersonic  regions  and  shock  waves,  a  large  number  of  grid  points  are  needed  to  capture 
the  flow  features.  This  is  especially  true  of  three  dimensional  transonic  flows.  Large 
computational  grids  require  supercomputers  with  large  memory  and  computing  speeds. 

Figure  2  shows  a  two-dimensional  cross-section  of  the  full  3D  grid  in  the  vicinity  of 
the  projectile.  The  surface  points  on  the  projectile  els  well  as  the  longitudinal  grid  point 
distribution  in  the  wake  are  selected  using  an  interactive  design  program.  The  grid  for  the 
body  region  of  the  projectile  was  computed  using  a  hyperbolic  grid  generation  program 
(Nietubicz,  Heavey,  Etnd  Steger  1982).  Each  grid  section  was  obtained  separately  and  then 
appended  to  provide  the  full  grid.  The  full  grid  is  split  into  three  zones,  one  in  front  of 
the  nose,  one  on  the  projectile  and  the  third  one  in  the  wake  or  base  region.  These  grids 
consist  of  30x79,  175x50  and  59x99  grid  points,  respectively.  The  3D  grids  consist  of  39 
points  in  the  circumferential  direction  for  the  symmetric  (non-spinning)  case  and  72  points 
for  spinning  case  with  no  symmetry.  An  expanded  view  of  the  wake  region  grid  is  shown 
in  Figure  3.  It  shows  the  longitudinal  grid  clustering  on  the  boattail  and  near  the  base 
corner  of  the  projectile.  The  grid  points  axe  also  clustered  near  the  projectile  surface  to 
capture  the  viscous  effects  in  the  boundary  layer.  For  simplicity,  this  clustering  has  been 
extended  from  the  base  corner  downstream  into  the  wake  region. 

The  current  problem  of  interest  is  the  effect  of  spin  on  the  projectile  aerodynamics. 
The  large  grids  required  for  accurate  numerical  simulations  of  three  dimensional  transonic 
flows  necessitates  the  use  of  a  Cray  class  supercomputer.  All  numerical  calculations  for  the 
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present  study  were  made  on  the  Cray-2  supercomputer  at  the  Ballistic  Research  Laboratory 
(BRL).  The  numerical  simulation  of  the  non-spinning  (symmetric)  cases  requires  about  44 
million  words  of  central  core  memory  on  the  Cray-2  whereas  full  3D  calculations  for  the 
spinning  cases  with  no  plane  of  symmetry  requires  about  82  million  words  of  memory. 
Each  of  the  numerical  calculations  for  the  symmetric  case  takes  about  30  to  40  hours  of 
computer  time.  The  time  required  to  get  the  final  result  for  the  spinning  case  is  about  two 
times  that  of  the  symmetric  case. 


Figure  2.  Computational  grid  (wind  and  lee  side). 


5.  RESULTS 

Numerical  computations  have  been  made  for  both  spinning  and  non-spinning  cases.  All 
the  computations  have  been  run  at  Mo0=  0.94  and  a  =  0,4  and  10°.  Atmospheric  flight 
conditions  were  used.  The  three  plane  version  of  the  3D  code  was  run  for  the  zero  degree 
angle  of  attack  case.  For  the  non-spinning  case,  39  planes  were  used  in  the  circumferential 
direction.  Two  end  planes  were  used  to  specify  symmetric  boundary  conditions  with  the 
remaining  37  used  for  the  circumferential  half  plane.  The  number  of  planes  used  for  the 
full  3D  spinning  case  was  72.  For  all  these  grids  the  circumferential  grid  spacing  was  5°. 
For  the  zero  angle  of  attack  case,  the  grid  was  rotated  circumferentially  5°  on  either  side 
of  the  mid  plane.  This  provided  the  3  planes  needed  in  the  code  to  use  centred  finite 
differences  in  the  circumferential  direction.  In  each  case,  the  solution  was  marched  from 
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free  stream  conditions  everywhere  until  the  final  converged  solution  was  obtained.  The 
results  are  now  presented  in  the  form  of  surface  pressure  coefficients,  and  aerodynamic 
force  and  moment  coefficients.  Comparison  of  the  computed  results  are  made  with  the 
available  experimental  data  (Miller  and  Molnar  1986).  The  computed  force  and  moment 
coefficients  are  also  compared  with  another  set  of  experimental  data  (Platou  and  Nielsen 
1973)  where  forces  and  moments  were  obtained  using  a  force  and  moment  balance. 

A  few  qualitative  results  axe  presented  next.  Figures  4,  5,  and  6  show  typical  Mach 
number  contour  plots  for  Moo=0.94  at  a  =  0,4,  and  10°  for  the  wind-side  and  the  lee- 
side.  Mixed  regions  of  subsonic  and  supersonic  flow  have  been  determined.  The  subsonic 
regions  are  found  in  front  of  the  flat  nose,  in  the  boundary  layer  and  in  the  near  wake. 
The  flow  regions  near  the  ogive-cylinder  junction  as  well  as  on  the  boattail  are  supersonic. 
Other  features  to  observe  are  flow  expansion  at  the  ogive-cylinder  and  cylinder-boattail 
corners  which  are  followed  by  shock  waves  (coalescence  of  contour  lines)  on  the  cylinder 
and  boattail  sections.  For  the  zero  degree  angle  of  attack  case  (see  Figure  4),  the  flow 
field  was  computed  only  for  the  top  half  plane  and  the  flow  field  has  been  mirrored  to 
emphasize  the  symmetric  location  of  the  shock  waves  along  the  body.  As  the  angle  of 
attack  is  increased  to  4  degrees  (Figure  5)  the  flow  field  becomes  asymmetric.  The  shock 
waves  are  asymmetrically  located  (the  wind-side  shocks  being  closer  to  the  base  than  the 
lee-side).  This  asymmetry  is  more  pronounced  with  a  further  increase  in  angle  of  attack  to 
10  degrees.  Figure  6  clearly  shows  the  asymmetrically  located  shock  waves  for  this  case. 


Figure  7  shows  the  longitudinal  surface  pressure  distribution  for  the  zero  degree  angle 
of  attack  case.  Figure  8  and  Figure  9  show  the  longitudinal  surface  pressure  distribution 
for  both  the  wind-side  and  lee-side  at  a  =  4.0°  and  a  =  10.0°,  respectively.  In  all  these 
cases  the  computed  surface  pressures  are  found  to  be  generally  in  good  agreement  with 
the  experimental  data  (Miller  and  Molnar  1986).  The  flow  expansions  at  the  ogive  and 
boattail  corners  axe  indicated  by  the  pressure  drop  at  these  junctions.  This  is  followed  by  a 
sharp  increase  in  pressure  at  the  shocks  (on  the  cylinder  and  boattail).  On  the  nose  section 
of  the  projectile,  for  the  angle  of  attack  cases,  the  pressure  on  the  wind-side  (  <f>  —  0.0°  ) 
is  greater  than  the  pressure  acting  on  the  lee-side  (  <f>  =  180.0°  ).  The  reverse  is  true  on 
the  boattail  section.  As  shown  in  Figure  8  for  a  =  4.0°,  the  change  in  pressure  between 
the  wind-side  and  the  lee-side  on  the  cylindrical  section  is  rather  small.  As  angle  of  attack 
is  increased  to  a  =  10.0°,  the  difference  in  the  wind-side  and  the  lee-side  pressures  gets 
larger  especially  on  the  ogive  and  boattail  sections  (  Figure  9  ). 


Figure  4.  Mach  number  contours,  M(X>= 0.94,  a  =  0.0°. 

The  next  series  of  figures  show  the  circumferential  pressure  distribution  for  both  non¬ 
spinning  and  spinning  cases.  The  results  for  the  non-spinning  cases  are  shown  first.  Here, 
<t>  =  0.0°  corresponds  to  the  wind-side  while  the  lee-side  is  at  <f>  =  180.0°.  Although, 
the  computations  are  performed  only  for  the  half  body,  the  computed  circumferential 
pressure  distributions  have  been  mirrored  and  are  shown  for  the  full  body  i.e.  for  <j>  =  0  to 
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Figure  7.  Longitudinal  surface  pressure  distribution,  Mqq— 0.94,  a  —  0.0° 
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Longitudinal  surface  pressure  distribution,  Moo=0.94,  a  =  4.0 
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360.0°.  Figure  10  shows  the  surface  pressure  distribution  for  a  longitudinal  position  X/D 
=  4.76  which  is  on  the  cylindrical  portion  of  the  projectile  ahead  of  the  boattail  corner. 
Figures  11  through  14  show  the  surface  pressure  distribution  at  X/D  =  5.11,  5.30,  5.43, 
and  5.50,  respectively.  These  longitudinal  positions  axe  all  located  on  the  boattail  of  the 
projectile.  These  figures  show  the  effect  of  angle  of  attack  on  the  pressure  distribution. 
The  computed  pressures  are  compared  with  experimental  data  (Miller  and  Molnar  1986) 
in  each  of  these  figures.  As  seen  in  Figure  10  for  X/D  =  4.76,  the  variation  of  pressure 
with  4>  is  almost  negligible  for  a  =  4.0°.  However,  at  a  =  10.0°  the  pressure  first  drops  and 
then  increases  in  going  from  the  wind-side  to  the  lee-side.  The  computed  pressures  are  in 
very  good  agreement  with  the  data.  For  X/D  =  5.11  (Figure  11),  there  exists  a  noticeable 
circumferential  pressure  variation  for  both  a  =  4.0°  and  a  =  10.0°.  The  extent  of  the 
circumferential  variation  in  pressure  increases  with  locations  further  aft  on  the  boattail  as 
shown  in  Figure  11  through  Figure  14.  For  a  =  4.0°,  the  pressure  variation  is  small  near 
the  wind-side  and  then  the  pressure  rises  as  the  lee-side  is  approached.  For  a  =  10.0°, 
however,  the  pressure  first  drops  and  then  increases  sharply  as  the  lee-side  is  approached. 
Computed  pressure  peaks  occur  on  the  lee-side  for  both  cases.  The  computed  pressures 
are  generally  in  good  agreement  with  the  experimental  data.  Also,  as  seen  in  these  figures, 
large  variations  in  surface  pressures  occur  with  increase  in  angle  of  attack. 


<p  (degrees) 

Figure  10-  Circumferential  surface  pressure  distribution,  P  =  0,  M00= 0.94, 
X/D=4.76. 
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Figure  11.  Circumferential  surface  pressure  distribution,  P  =  0,  MCO=0.94, 

X/D=5.11. 
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Figure  12.  Circumferential  surface  pressure  distribution,  P  =  0,  Mqq 
X/D=5.30. 


Figure  13.  Circumferential  surface  pressure  distribution,  P  =  0,  M ^ 
X/D=5.43. 
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Figure  14.  Circumferential  surface  pressure  distribution,  P  —  0,  Ma O=0.94. 

X/D=5.50. 

The  results  for  the  spinning  cases  are  shown  next.  Here  again  the  circumferential 
pressure  distributions  are  shown  at  the  same  longitudinal  positions  X/D  =  4.76.  5.11. 
5.30,  5.43,  and  5.50.  Figure  15  shows  the  pressure  distribution  at  X/D  =  4.76.  As  angle 
of  attack  is  increased,  variation  of  pressure  with  4>  starts  showing  up.  The  agreement  of 
the  computed  pressure  with  the  data  is  good  for  a  =  10.0°  .  Some  discrepancy  exists  for 
the  a  =  4.0°  case.  The  computed  result  for  this  low  angle  of  attack  case  has  been  checked 
and  should  be  at  least  as  accurate  as  the  high  angle  of  attack  case  at  X/D  =  4.76  (  on 
the  cylinder).  The  discrepancy  at  this  particular  location  is  not  a  result  of  computational 
inaccuracy.  Possibly,  the  data  at  this  location  is  suspect.  Figure  16  through  Figure  19 
show  the  surface  pressure  distributions  on  the  boattail  at  X/D  =  5.11,  5.30,  5.43,  and  5.50. 
respectively.  These  figures  again  show  the  large  variation  of  pressure  with  increase  in  angle 
of  attack  and  with  distance  downstream  on  the  boattail.  The  agreement  of  the  computed 
pressures  with  the  experimental  data  is  not  as  good  as  it  was  found  for  the  non-spinning 
case. 

Figures  20  and  21  show  the  effect  of  spin  on  the  pressure  distribution  at  X/D=4.76 
for  a  =  4.0°  and  a  =  10.0°,  respectively.  As  expected,  the  spin  induces  small  changes  in 
the  pressure  distribution.  For  the  ten  degree  angle  of  attack  case,  the  pressure  is  reduced 
slightly  with  spin.  This  is  predicted  by  the  computations  and  is  also  observed  in  the 
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Figure  15.  Circumferential  surface  pressure  distribution,  P  =  4900  rpm,  M^=0.94, 
X/D=4.76. 
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Figure  16.  Circumferential  surface  pressure  distribution,  P  =  4900  rpm.  .A/x-  =0.94. 
X/D=5.11. 
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Figure  17.  Circumferential  surface  pressure  distribution,  P  =  4900  rpm,  A/^,=0.94. 
X/D=5.30- 
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Figure  18.  Circumferential  surface  pressure  distribution.  P  =  4900  rpm.  M, x>=0.94, 
X/D=5.43. 
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Figure  19.  Circumferential  surface  pressure  distribution,  P  =  4900  rpm,  Mp^O.94, 
X/D=5.50. 

experimental  data  (Miller  and  Molnar  1986).  For  the  four  degree  angle  of  attack  case, 
the  computation  again  shows  the  same  trend;  however,  the  data  show  the  reverse  trend 
from  4>  =  40.0°  to  <f>  =  160.0°.  Figures  22  and  23  show  the  effect  of  spin  on  the  pressure 
distribution  at  X/D=5.43  for  a  =  4.0°  and  a  =  10.0°,  respectively.  For  the  zero  spin 
cases  shown  here,  the  pressure  distribution  is  symmetric  about  4>  =  180.0°.  As  shown  in 
Figure  22  for  a  —  4.0°,  the  agreement  of  computed  pressure  with  data  for  the  spinning 
case  is  good  except  for  <f>  =  240.0°  to  (j>  =  360.0°.  As  shown  in  Figure  23  for  the  a  =  10.0° 
case,  the  trend  of  computed  pressure  with  spin  seems  to  agree  fairly  well  with  the  data; 
however,  discrepancies  do  exists  in  the  quantitative  agreement.  For  this  high  angle  of 
attack  case,  the  experimental  pressure  distribution  also  shows  a  second  pressure  increase 
from  <p  =  140.0°  to  4>  =  170.0°  near  the  lee-side  which  is  not  predicted  in  the  computed 
pressures. 

The  entire  flowfield  over  the  projectile  including  the  base  region  is  computed.  There¬ 
fore,  the  computed  results  include  any  upstream  influence  the  base  region  flow  may  have 
on  the  boattail  flowfield.  Surface  pressures  including  the  base  pressure  and  the  viscous 
stresses  are  known  from  the  computed  flow  field  and  can  be  integrated  to  give  the  aerody¬ 
namic  forces  and  moments. 
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Figure  21.  Circumferential  surface  pressure  distribution,  MCO=0.94,  a  = 


X/D=4.76. 
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Figure  22.  Circumferential  surface  pressure  distribution,  M00= 0.94,  a  = 


X/D=5.43. 
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Figure  23.  Circumferential  surface  pressure  distribution,  =0.94,  a  = 


X/D=5.43. 
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Figure  24  and  Figure  25  show  the  effect  of  angle  of  attack  on  normal  force  and  pitching 
moment  coefficient,  respectively.  Computed  results  are  compared  with  the  experimental 
data  (Miller  and  Molnar  1986).  In  addition,  these  figures  include  Platou  data  (1973) 
which  were  obtained  using  a  force  and  moment  balance.  As  seen  in  both  the  computed 
and  experimental  data,  the  effect  of  spin  on  these  coefficients  is  small.  The  normal  force 
increases  and  the  pitching  moment  coefficient  decreases  with  angle  of  attack.  This  is 
the  correct  trend  seen  in  both  the  computed  results  and  data.  As  shown  in  Figure  24 
for  normal  force  coefficient,  discrepancy  of  the  order  of  10  —  15%  does  exist  between  the 
computed  result  and  experimentally  obtained  coefficient  data  at  a  =  4.0.  The  agreement 
is  much  better  at  a  =  10.0.  The  discrepancy,  particularly  for  the  non-spinning  cases  at 
low  angle  of  attack,  is  difficult  to  explain  since  the  computed  surface  pressures  are  in 
good  agreement  with  the  experimentally  measured  pressures  (Miller  and  Molnar  1986). 
Note  that  these  experimental  measurements  were  obtained  at  limited  positions  along  the 
body  (practically  no  points  on  the  ogive  and  only  five  points  on  the  boattail).  As  seen  in 
Figure  25,  the  pitching  moment  coefficients  obtained  experimentally  from  measured  surface 
pressures  (Miller  and  Molnar  1986)  are  smaller  than  Platou  data  (1973).  The  computed 
results  lie  in  between  the  two  sets  of  experimental  results. 
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Figure  24.  Normal  force  coefficient  vs.  angle  of  attack,  Mqq =0.94. 


An  aerodynamic  coefficient  which  is  of  primary  concern  is  the  Magnus  force  coefficient, 
Cy  •  Figure  26  shows  the  Magnus  (side)  force  as  it  develops  along  the  projectile  for  the 
spinning  case  at  a  =  4.0  and  10.0°.  For  both  of  these  cases,  the  side  force  contribution 
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Figure  25.  Pitching  moment  coefficient  vs.  angle  of  attack,  Moa= 0.94. 


due  to  viscous  stresses  is  negligible.  The  computed  Magnus  force  shown  in  this  figure  then 
comes  primarily  from  the  integration  of  surface  pressures.  The  computed  results  are  shown 
in  solid  and  dash  lines  .  The  open  symbols  represents  the  experimental  results  obtained 
from  measured  surface  pressures  (Miller  and  Molnax  1986).  The  closed  symbols  represent 
wind  tunnel  measurements  (Platou  and  Nielsen  1973)  where  Magnus  force  and  moment 
were  obtained  using  a  balance.  For  the  a  =  4.0°  case,  both  the  computed  result  as  well 
as  the  data  show  negligible  contribution  from  the  ogive  section  of  the  projectile  and  are 
in  good  agreement  with  each  other.  The  side  force  drops  right  after  the  ogive-cylinder 
junction  and  then  begins  to  level  out.  Again  both  results  show  this  and  the  agreement 
is  good  until  X/D  =  4.0,  half  way  through  the  cylindrical  section.  From  here  on  the 
computed  result  shows  a  gradual  drop  in  side  force  for  the  rest  of  the  cylindrical  section 
while  the  experimental  result  shows  a  shaxp  decline.  At  the  boattail  junction,  there  is  a 
small  increase  in  this  force.  This  is  followed  by  a  shaxp  drop  for  almost  all  of  the  boattail 
except  near  the  base  corner  where  it  increases  somewhat.  This  trend  can  also  be  observed 
in  the  computed  result  and  was  shown  in  the  earlier  computations  (Nietubicz,  Sturek,  and 
Heavey  1983).  The  final  value  of  the  side  force  agrees  well  with  the  experimental  results 
at  a  =  4.0°.  For  the  a  =  10.0°  case,  the  trend  of  the  side  force  distribution  observed  in 
the  experimental  result  (Miller  and  Molnax  1986)  can  also  be  seen  in  the  computed  result. 
However,  there  is  discrepancy  between  both  of  the  experimental  results  as  well  a is  between 
the  computation  and  experiment. 
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Figure  26.  Development  of  Magnus  force  along  the  projectile,  Moo=0.94. 


Another  aerodynamic  coefficient  which  is  very  important  from  dynamic  stability  point 
of  view  is  the  Magnus  (Yawing)  moment  coefficient,  Cn.  Figure  27  shows  the  Magnus  mo¬ 
ment  coefficient  for  the  projectile  as  a  function  of  angle  of  attack.  This  moment  coefficient 
is  referenced  to  the  center  of  gravity  of  the  projectile  which  is  located  at  3.506  calibers 
from  the  nose.  The  computed  result  is  shown  by  the  solid  line.  The  circles  represent  the 
experimental  results  (Miller  and  Molnar  1986)  and  the  boxes  represent  Platou  data  (1973) 
obtained  using  a  force  and  moment  balance.  In  general,  the  computed  result  is  in  good 
agreement  with  the  experimental  data  of  Miller  and  Molnar  (1986).  For  the  a  =  4.0°  case, 
both  the  computed  result  as  well  as  the  data  are  in  reasonable  agreement  with  each  other. 
However,  at  a  =  10.0  there  is  some  discrepancy  or  disagreement  even  between  the  two 
experimental  results. 

If  the  dependence  of  the  Magnus  moment  is  cubic  in  yaw  level  or  angle  of  attack,  the 
nonlinear  variation  of  the  moment  coefficient  is  of  the  general  form: 


=  Cnao  +  C2S2  (7) 

where  is  the  zero-yaw  value  of  Magnus  moment  coefficient,  C2  is  the  associated  cubic 
coefficient  and  6  =  sina.  This  analysis  of  the  nonlinear  Magnus  moment  yields  the  follow¬ 
ing  result. 
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Figure  27.  Magnus  moment  coefficient,  MQO= 0.94. 


Table  1.  Analysis  of  Nonlinear  Magnus  moment 


Coefficient 

Computation 

Miller’s 

experiment 

Platou’s 

data 

.196 

.210 

.245 

c2 

-.48 

-.55 

As  shown  in  Table  1,  the  computed  zero- yaw  Magnus  moment, C„Qo,  agrees  to  within  about 
6%  and  the  computed  cubic  coefficient  agrees  to  within  about  12%  of  Miller’s  (1986) 
experimental  data  for  the  exact  same  shape  and  flow  conditions.  It  should  be  noted  that 
this  result  represents  the  first  numerical  prediction  of  a  nonlinear  Magnus  moment  at  a 
transonic  velocity. 

As  discussed  earlier  for  this  high  angle  of  attack  (a  =  10.0°),  the  experimental  pressure 
distribution  shows  a  second  pressure  increase  (hump)  from  <f>  =  140.0°  to  (f>  =  170.0°  near 
the  lee-side  which  is  not  predicted  in  the  computed  pressures  (see  Figure  22).  It  is  well 
known  that  cross  flow  separation  occurs  on  the  projectile  at  high  angles  of  attack.  One  way 
to  possibly  improve  the  numerical  predictions  for  a  high  angle  of  attack  case  is  to  compute 
the  lee-side  vortical  flow  structure  accurately.  The  turbulence  level  predicted  by  Baldwin- 
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Lomax  algebraic  turbulence  model  can  be  large  in  this  vortical  flow  region.  Limiting  the 
outer  turbulence  length  scale  in  the  model  in  this  region  has  led  to  improved  numerical 
results  (Schiff,  Degani,  and  Cummings  1989)  for  supersonic  and  subsonic  flows.  This  idea 
has  been  used  here  and  a  solution  has  been  obtained  for  a  =  10.0°.  Figure  28  and  Figure 
29  show  the  circumferential  surface  pressure  distribution  at  two  longitudinal  positions,  one 
on  the  cylinder  (X/D  =  4.76)  and  the  other  on  the  boattail  (X/D  =  5.50).  These  figures 
show  the  result  obtained  by  the  modified  turbulence  model.  This  result  is  compared  with 
that  obtained  by  the  original  turbulence  model  as  well  as  the  experimental  data.  As 
shown  in  Figure  28  for  X/D  =  4.76,  the  result  obtained  with  the  modified  turbulence 
model  agrees  better  with  the  data  and  also  contains  the  pressure  hump  near  the  lee-side. 
At  X/D  =  5.50  (see  Figure  29)  the  agreement  of  the  computed  result  with  the  data  is 
not  as  good.  However,  the  same  trend  of  the  variation  of  pressure  with  <j)  observed  in  the 
experimental  data  (Miller  and  Molnar  1986)  has  been  clearly  predicted  in  the  computation 
with  the  modified  turbulence  model.  This  includes  all  the  pressure  variations  near  the  lee- 
side.  The  aerodynamic  forces  and  moments  (including  the  Magnus  force)  obtained  from 
these  pressures  were  found  to  change  very  little  from  the  results  obtained  with  the  original 
turbulence  model. 


cp  (degrees) 

Figure  28.  Circumferential  surface  pressure  distribution,  M00= 0.94,  a  =  10.0°, 
X/D=4.76. 
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Figure  29.  Circumferential  surface  pressure  distribution,  M00= 0.94,  a  =  10.0°, 
X/D=5.50. 

6.  CONCLUDING  REMARKS 

A  zonal,  implicit,  time-maxching  Navier-Stokes  computational  technique  has  been  used 
to  compute  three  dimensional  transonic  flowfields  over  a  secant-ogive-cylinder- boattail  pro¬ 
jectile.  Flowfield  computations  have  been  performed  at  M  =  0.94  for  spin  rates  of  0  and 
4900  rpm  and  at  angles  of  attack,  a  =  0, 4,  and  10°  . 

Numerical  results  show  the  details  of  the  flowfield  such  as  Mach  number  contours  and 
surface  pressure  distributions.  Computed  surface  pressures  are  compared  with  available 
experimental  data  for  the  same  conditions  and  the  same  configuration.  Computed  re¬ 
sults  show  the  large  circumferential  pressure  variations  over  the  boattail  region  as  well  as 
the  non-linear  effect  of  the  angle  of  attack.  Aerodynamic  force  and  moment  coefficients 
(normal  force,  pitching  moment  coefficient,  Magnus  force  and  Magnus  moment)  have  been 
obtained  from  the  computed  pressures  and  have  been  compared  with  the  data.  Reasonable 
agreement  has  been  obtained  for  the  a  =  4°  case  for  all  these  coefficients.  The  comparison 
is  less  satisfactory  at  a  =  10°  particularly  for  the  Magnus  force.  A  modified  turbulence 
model  was  used  for  this  case  which  improved  the  pressure  comparison  somewhat.  However, 
it  did  not  change  the  force  and  moment  results  significantly.  It  should  also  be  noted  that 
for  the  high  angle  of  attack  case  both  experimental  results  do  not  agree  with  each  other. 
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Analysis  of  non-linear  Magnus  moment  was  made  to  obtain  both  the  zero-yaw  Magnus 
moment  as  well  as  the  cubic  coefficient.  The  computed  zero-yaw  Magnus  moment  agrees 
to  within  about  6%  and  the  computed  cubic  coefficient  agrees  to  within  about  12%  of 
Miller’s  (1986)  experimental  data  for  the  exact  same  shape  and  flow  conditions.  It  should 
be  noted  that  this  result  represents  the  first  numerical  prediction  of  nonlinear  Magnus 
moment  at  a  transonic  velocity. 

The  accuracy  in  the  experimentally  obtained  forces  and  moments  for  the  spinning  case  is 
largely  unknown.  Also,  it  is  not  clear  from  the  experiments  whether  the  boundary  layer  was 
fully  turbulent  or  if  transition  occurred  on  the  projectile.  The  boundary  layer  was  assumed 
to  be  turbulent  in  the  computations.  The  agreement  or  disagreement  of  the  computed 
forces  and  moments  with  the  experimental  results  must  be  viewed  in  that  light.  The  side 
force  is  very  small  and  it  is  very  difficult  to  obtain  experimentally  or  computationally. 
The  type  of  agreement  achieved  between  the  computation  and  experiments,  par*rc”1arly 
at  low  angle  of  attack  is  very  encouraging.  Additional  code  validation  is  required  to 
better  determine  the  accuracy  of  the  predicted  side  force  and  moment  by  comparison 
with  free  flight  results.  This  would  involve  full  3D  computations  at  a  low  angle  of  attack 
for  atmospheric  flight  conditions.  Future  efforts  will  then  include  running  calculations  at 
various  transonic  Mach  numbers  to  determine  the  critical  aerodynamic  behavior  in  Magnus 
force  and  moment. 
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speed  of  sound 

specific  heat  at  constant  pressure 
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Magnus  or  side  force  coefficient 

normal  force  coefficient 

pitching  moment  coefficient 

projectile  diameter 

total  energy  per  unit  volume 
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Prandtl  number 
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vector  of  dependent  variables 
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vector  containing  viscous  terms 
time 

temperature 

axial,  circumferential,  and  normal  velocity  components  of  the  Navier-Stokes  equations 
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a  Angle  of  attack 

7  ratio  of  specific  heats 

k  molecular  and  turbulent  thermal  conductivity 
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p  density 

<t>  circumferential  angle 
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oo  free  stream  conditions 
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